Third edge for a graphene nanoribbon: A tight-binding model calculation. 
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The electronic and transport properties of an extended linear defect embedded in a zigzag nanorib- 
bon of realistic width are studied, within a tight binding model approach. Our results suggest that 
such defect profoundly modify the properties of the nanoribbon, introducing new conductance quan- 
tization values and modifying the conductance quantization thresholds. The linear defect along the 
nanoribbon behaves as an effective third edge of the system, which shows a metallic behavior, giving 
rise to new conduction pathways that could be used in nanoscale circuitry as a quantum wire. 
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I. INTRODUCTION 

Graphene, a unique purely two dimensional system, 
evolved already into a landmark in condensed mater 
physics, since its first successful isolation [T] and the ini- 
tial establishment of their related electronic properties [2]. 
More recently, it has been addressed as to initiate a new 
generation of electronic devices [3], since the obtention 
of graphene nanoribbons 4 and the cutting out of more 
complex structures by means of atomic force microscopy 
nanolitography [5J ^.Nevertheless, many applications re- 
quire control at atomic scale [7j , which means that a con- 
flicting scenario has to be still overcome. On one side, for 
instance, clear single electron charging has been observed 
in graphene quantum dots [5] , but on the other side, there 
seem to be limitations to the use of graphene as ballis- 
tic nanowires due to still unavoidable disorder effects [9J . 
In this context, the introduction of localized defects has 
been suggested as one accessible way to modify the elec- 
tronic properties of graphene [TUHUI- Alternatively, the 
perspective of inducing extended defects, in particular 
linear defects, that could even assemble involved crossed 
structures over a large area on graphene, started to be 
discussed as a candidate for monolithic circuitry on a 
nanoscale [T31 [H] ■ Besides first principle studies on the 
stability of such structures and the associated density of 
states, systematic evaluation of their related transport 
properties have started to be undertaken [15] . 

A promising achievement is the observation of an ex- 
tended linear defect [16] on a graphene layer grown on 
a metallic substrate. This defect is created by alternat- 
ing Stone-Thrower- Wales (STW) defect Q3 [T5] and di- 
vacancies, leading to a pattern of repeating paired pen- 
tagons and octagons, as shown in Fig 1. It has been 
shown, combining first principle calculations and experi- 
mental scanning probe microscopy observations, that this 
defect acts as a one dimensional metallic wire JT3J [TB] , 
widening a road towards a graphene based electronics. 

Here we focus on a tight binding model calculation 
to scrutinize the conductance of such a linear defect em- 
broidered along a graphene nanoribbon with outer zigzag 




FIG. 1: Geometry of the extended defect considered. 



edges. After validating the model by ensuring that the 
main features of the electronic properties of the linear 
defect are captured by a nearest neighbor tight binding 
Hamiltonian, the conductance is obtained within a Lan- 
dauer Buttiker approach. The complex structure of the 
defect leads to sharp anti-resonances in the lower conduc- 
tance plateaus, which, on the other hand, may either oc- 
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cur at even multiples of Go = — , rather unexpected for 
graphene systems or may define new conductance plateau 
thresholds, acting effectively as a third edge in a zigzag 
nanoribbon. The present results are relevant for the ex- 
perimental context in which the controlled fabrication of 
extended defects [18] may be conbined with the nanoli- 
tographyc definition of nanoribbons [6] . 

In section II we briefly sketch our model approach and 
proceed by validating it in respect to previous first prin- 
ciple calculations [131 [TB] . In section III the electronic 
band structure of an infinite zizgag nanoribbon with and 
embedded linear defect is shown. In section IV the con- 
ductance , as well as the real space current density distri- 
bution modification due to the presence of a linear defect 
is discussed, while in section V the final comments and 
experimental perpsectives are placed. 



FIG. 2: Density of states (DOS) for an extended defect in an 
infinite graphene sheet. 



II. HEURISTIC MODEL FOR AN EXTENDED 
DEFECT 

We model the electronic structure by a nearest neigh- 
bor tight-binding Hamiltonian for an hexagonal lattice: 

H = t (cUj+cta) (1) 

<i,j> 

where Cj is the fermionic operator on site i. We model the 
planar arrangement of the atoms in the extended linear 
defect as shown in Fig. 1 by rearranging which atoms 
are connected by nearest-neighbor hopping parameters 
(i). This type of defect shows atomic separations rang- 
ing from 1.38 A to 1.44 A and conserves the coordination 
number |16j . which implies variations lesser than 5% in 
t, suggesting that t values could be considered nearly 
unaffected in respect to a defect free graphene nanorib- 
bon (for which t « 2.7 eV). In order to compare the 
results from our model to previous calculations O US] , 
we calculate the density of states (DOS) for an infinite 
graphene sheet (18200 atoms in the unit cell) with the 
extended defect. The result is shown in Fig. 2 where a 
peak around the Dirac point (E/t =0) is developed, in- 
dicating the metallic character and the one dimensional 
nature of the defect at this energy, in accordance to the 
already mentioned first principles calculations [13] I16j . 
The conductance G(E) is evaluated within the Landauer- 
Biittiker formalism, G(E) = G T(E), where G = ^ is 
the conductance quantum and T(E) is the transmission 
function between the contacts and is evaluated by 20 : 

Tpq = Tr[T p G pq T q G pq ] (2) 

where G r pq is the Green function between the contact 
p and q evaluated thorough the recursive lattice Green's 
function technique 21, 22 . The contact broadening func- 
tion [23] r p ( g ) = i(S p (g) - where E p(g ) is the self- 
energy of the contact, arises from the interaction of the 
central region with the semi infinite contacts. In order 
to calculate this term, it is necessary to know the Green 
function of the contact, also obtained numerically [2~4] . 
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FIG. 3: (Color online) (a) Geometry of zigzag graphene 
nanoribbon of width W with an extended defect at a dis- 
tance rf e from the upper edge, (b) Band structure of a zigzag 
nanoribbon of width W = 15.9 nm without any defect, (c) 
Band structure of a zigzag nanoribbon of W = 15.9 nm with 
an extended defect at d e = 8 nm .(d) Band structure of a 
zigzag nanoribbon of W = 15.9 with an extended defect at 
d e = 3.8 nm 



III. A WIRE EMBEDDED IN A ZIGZAG 
NANORIBBON 

The band structure of an infinite linear defect em- 
bedded along a zigzag graphene nanoribbon is obtained 
within the scheme illustrated in Fig. 3(a), where the 
unit cell used for the calculation is also shown. The ob- 
jectives here are twofold: on one hand we have to com- 
pare the band structure obtained within our model and 
for a nanoribbon with the band structure of a bulk infi- 
nite graphene sheet, as presented previously [16]. On the 
other hand, we want to capture the changes in the band 
structure induced by the linear defect looking forward 
the interpretation of the conductance features of a scat- 
tering region of finite length (where the defect is embed- 
ded) terminated by semi infinite zigzag nanoribbon con- 
tacts of the same width. The band structure of a zigzag 
nanoribbon W = 15.9 nm wide, without the extended de- 
fect, is depicted in Fig. 3(b): a flat band at zero energy, 
corresponding to edge states, is developed in the range 
— 27r/3 < k < 2tt/3. Notice that usually these flat bands 
are seen within 27r/3 < |fc| < 7r [7]. Our slightly different 
representation is due to a folding introduced by the fact 
that we use a unit cell twice as large as necessary for a 
defect less nanoribbon, in order to properly describe the 
defect afterwards. The presence of the extended defect 
in the nanoribbon does not alter the metallic character 
or the appearance of the edge states in the nanoribbon, 
but clearly (i) breaks the particle-hole symmetry |25j in 
a new context, due to the rebounding along the extended 
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defect and (ii) creates a new flat band at zero energy (now 
in a range around k — 0, in accordance to what has been 
obtained for the bulk case [16]) which bends upwards in 
energy as shown in Fig. 3(c), for an extended defect 
located in the middle (d e — 8 nm) of the nanoribbon. 
Mapping the charge density (not shown here) for differ- 
ent values of k related to this band, one could see that the 
wave function spreads out mainly along the edges with a 
significant contribution at the defect sites, a contribution 
that rapidly decays away from the defect region. This de- 
cay of the wave function away from the defect has been 
already measured [TCj and our findings within the present 
heuristic model are qualitatively in accordance with the 
experimental data. The degeneracy of the states at zero 
energy is broken at different k values depending on the 
width of the nanoribbon or by simply putting the defect 
closer to one edge, the case shown in Fig. 3(d), depicting 
the band structure when the extended defect is located 
at d e — 3.8 nm, i.e., roughly halfway towards an edge 
from the middle, in the same nanoribbon of W = 15.9 
nm. Comparing Figs. 3(c) and 3(d), the stronger warp- 
ing of the defect-like band is stronger for a defect line 
closer to the edge, but the effect seems to be qualita- 
tively the same, irrespective to the distance to one of the 
outer edges. More importantly is that, comparing the de- 
fect free nanoribbon, Fig. 3(b), with Fig. 3(c) and 3(d), 
one sees the clustering in pairs of the bulk bands, a fin- 
gerprint of a division of the original nanoribbon into an 
effective double ribbon in parallel. These results for an 
infinite system constitute key ingredients to understand 
the new conductance features. 



IV. CONDUCTANCE IN THE PRESENCE OF 
THE EXTENDED DEFECT 

A. Modified conductance quantization 

The conductance is calculated for nanorribons of dif- 
ferent lengths L, containing longitudinal linear defects. 
Recalling, for sake of clearness, these systems are at- 
tached to two semi-infinite contacts of ideal graphene 
zigzag nanoribbons of the same width 15.9 nm. The con- 
ductance of a perfect nanoribbon shows quantized con- 
ductance plateaus following an odd sequence specific rule 
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G(E) = G (2n + l) 



(3) 



while the conductance plateaus onset energies are given 
by |2E]: 



E n . , ,1 , V^air 



(4) 



where Go = n indicates the number of the con- 
ducting channel, a = 2.4 A is the graphene lattice con- 
stant and W is the nanoribbon width. For W = 15.9 
nm & R3 ±0.06 and & f=s ±0.11. The ideal conductance 
line shapes for defect free nanoribbon for the present pa- 
rameters are plotted in Fig. 4 [(a)-(d)] as dashed lines. 
However, in Figs. 4(e) and (f), the dotted lines are for the 
conductance of a nanoribbon with an infinite linear de- 
fect, corresponding to the band structure shown in Fig. 3 
(c) and ((d), respectively. The actual conductance, i.e., 
including the finite longitudinal defects, is depicted as 
continuous lines. In what follows we discuss the impor- 
tant and also departures from the transport of an ideal 
zigzag nanoribbon. First of all, Fig. 4 shows nanoribbons 
of different length (with the linear length of the order of 
the nanoribbon), namely L = 50 nm and L = 375 nm, 
and different defect to edge distances: d e — 8 nm, i.e. in 
the middle of the ribbon, and d e = 3.8 nm, half way from 
the middle towards one of the edges. 

Let us consider first Fig. 4(a), a shorter defect in the 
middle of the ribbon, as a starting point. Initially, one 
has to recall the band structures in Fig. 3 and the particle- 
hole symmetry breaking, with one of fiat bands bend- 
ing downwards as in ideal ribbons, while a defect related 
band bends upwards together with a edge related band. 
Therefore the conductance shows two different kinds of 
oscillations: there are no qualitative changes in the first 
hole-like plateau region, the oscillations shown in this 
energy range are simply Fabry-Perot interference effects 
due to the broken translational symmetry introduced by 
the defect. On the other hand, the first electron-like 
plateau shows not interference oscillations but strong 
anti-resonances [33], due to the coupling of extended 
states at the edges with localized ones in the defects, as 
discussed elsewhere for a double vacancy system j3U] . No- 
tice that these anti resonances change into Fabry-Perot 
oscillations in the second electron-like plateau. This is 
associated to a change from localized to extended char- 
acter of the defect related states in this energy range, 
where, besides the interference oscillations, the conduc- 
tance maxima reaches 2Go instead of the expected 3Go 
for clear graphene nanoribbons. 

Looking at a longer defect system, Fig. 4(b), one 
sees a clearer definition of the even electron conductance 
plateau, besides the increase in the number of interfer- 
ence oscillations, as well as anti resonances in all plateaus, 
simply due to increasing the length of the extended de- 
fect. Here one should pay attention to the third electron 
conductance plateau that evolves in a non-trivial way 
towards an expected 5Go value. Besides the strong oscil- 
lations, analogous to the already discussed ones in lower 
plateaus, it can be seen that this plateau splits into two 
sub plateaus, one reaching 3G and only the second ap- 
proaching 5Go- This plateau splitting, occurring also at 
other higher plateaus and in the hole side of the spec- 
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FIG. 4: (Color Online) Dashed lines correspond to the con- 
ductance for a perfect zigzag graphene nanoribbon of width 
W = 15.9 nm. Solid lines correspond to the conductance for 
an extended defect inserted in a ribbon with (a) ribbon width 
W = 15.9 nm, defect length L = 50 nm and distance from 
defect to the upper edge d e = 8 nm. (b) W = 15.9 nm, 
L = 375 nm and d e = 8 nm. (c) W = 15.9 nm, L = 50 
nm and d e = 3.8 nm. (d) W = 15.9 nm, L = 375 nm and 
d e = 3.8 nm. (e) Dotted line correspond to the conductance 
of a zigzag nanoribbon, W = 15.9, with an infinite linear de- 
fect at d e — 8 nm. Continuous line W = 15.9 nm, defect 
length L — 50 nm and d e = 8 nm. (f) Dotted line corre- 
spond to the conductance of a zigzag nanoribbon, W = 15.9, 
with an infinite linear defect at d e — 3.8 nm. Continuous line 
W — 15.9 nm, defect length L = 50 nm and d e = 3.8 nm. 
Arrow refers to Fig. 5 



trum, is due to the band clustering in pairs, already dis- 
cussed in respect to Fig. 3. Indeed, the onsets of these 
sub plateaus are given by the onset of subsequent bands 
in Fig. 3. 

Moving now the defect towards one of the edges, new 
features appear, as observed in Figs. 4(c) and 4(d). The 
most striking one is the fact that the onset of second 
plateau onset, either for holes and electrons can not be 
associated to the corresponding ideal nanoribbon, as hap- 
pens for the linear defect in the middle. It is worth not- 
ing that these effective conductance plateau onsets cor- 
respond to a thinner nanoribbon, namely W = 12.1 nm 
wide by fitting the numerical results in Fig. 4 by eq. 



(4). This corresponds exactly to the distance from the 
linear defect to the more distant edge: 15.9 nm - 3.8 
nm = 12.1 nm, recalling that the width of the nanorib- 
bon is W — 15.9 nm and d e — 3.8 nm is the distance 
of the defect to one of the edges, for the results shown 
in Fig. 4(c) and Fig. 4(d). This behavior, potentially 
observable, is explained by the fact that for d e = 3.8 
nm, the linear defect is strongly coupled to the nearest 
edge, defining, therefore a thinner effective nanoribbon. 
Another striking effect is the wide anti resonance shown 
in the first hole-like plateau, which can also be directly 
associated to the band structure for the corresponding in- 
finite system, Fig. 3(d): the edge state band that bends 
downwards presents a second flat step at the energy of 
the anti-resonance. 

To further interpret our preceding findings about the 
conductance features due to the extended defect, a di- 
rect comparison to the conductanceof ribbons with an 
infinite linear defect, i.e., corresponding directly to the 
band structure in Fig. 3, becomes necessary. These results 
are given as dotted lines in Figs. 4(e) and 4(f), consider- 
ing both d e — 8 nm d e = 3.8 nm, respectively. Although 
some features, like the existence of a 2Gq electron-like 
plateau, are revealed, the finite length defect nanoribbon 
is indeed an involved system, that can not be completely 
understood by means of infinite counterparts. In partic- 
ular, although the prediction of a 2Go plateau already at 
the Dirac point, the actual conductance for low energy 
electrons reaches only Go, as expected from a defect free 
nanoribbon, actually the description of the leads con- 
nected to the defect region. Hence, the finiteness of the 
linear defects introduces unexpectancies. The supression 
of the conductance at the first electron-like plateau in 
Fig. 4(e) and 4(f) are due to backscattering of a defect 
related mode at the end of the connection to the leads. 



B. Metallic third edge 

To further investigate the linear defect characteristics, 
the real space current density distribution is depicted in 
units of Go, following the procedure sketched in [231131"]. 
Fig. 5. The arrow in Fig. 4(e) points to the first conduc- 
tance peak reaching 3Go at an energy (E/\t\ — 0.0651), 
associated to the real space current density distribution 
shown in Fig. 5(b). A high current density around the 
defect reveals its metallic character, although the current 
spreads out considerably across the ribbon, differently 
from defect free system, where the low energy current 
density is located mainly at the center of the nanoribbon 
[3"2"] . A closer look shows that the defect acts as a third 
edge, since its presence creates a higher current density in 
the middle of the two newly defined nanoribbons around 
the linear defect. Interference fringes on both sides of 
this central edge, are also appreciated, as well as in the 
region in contact with the pristine graphene nanoribbon 




FIG. 5: (Color Online) (a) Conductance of a zigzag ribbon 
width W = 15.9 nm with a linear defect of length L — 50 nm 
and distance from defect to the upper edge d e — 5.8 nm. (b) 
Real space current density distribution, in units of Go = , 
for the conductance peak indicated by the arrow in Fig. 4(e). 
(c) Real space current density distribution, in units of Go, for 
the conductance peak indicated by the arrow in Fig. 5(a). 

leads. A similar behavior is found (high current density 
located along the embeded linear structure) if the cur- 
rent density is plotted for the others conductance peaks 
in this electron plateau (not shown here). In Fig 5(a) 
the conductance of a nanoribbon of width W — 15.9 nm, 
d e = 5.8 nm and L = 50 nm, i.e., an off center embedded 
edge, is shown, where again the second electron plateau 
reaches the value of 2Gq. The arrow points to the first 
conductance peak defined at an energy of E/\t\ = 0.0662 
and the real space current density distribution at this en- 
ergy is shown in Fig 5(c). The high current density along 
the defect is appreciated as well as its role as a third edge, 
defining two clearly different regions with different cur- 
rent densities, loosing now the symmetry of the current 
density between the upper ribbon and the lower ribbon. 

One could envisage the defect as part of nanoscale cir- 
cuitry, such system can be created by scanning probe 
microscopy patterning a nanoribbon 5 in an appropri- 
ate region of a graphene sheet containing such a linear 
structure. If the linear defect is a rather ordered struc- 
ture |16j . the patterned outer edges will become quite 
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FIG. 6: (Color Online) (a) Average Conductance of a zigzag 
ribbon of width W = 15.9 nm, defect length L = 50 nm and 
distance from defect to the upper edge d e — 8 nm with edge 
disorder. Dashed line without edge disorder, (b) Average 
Conductance of a zigzag ribbon of width W = 15.9 nm, defect 
length L = 50 nm and distance from defect to the upper edge 
d e = 5.8 nm with edge disorder. Dashed line without edge 
disorder. 

disordered. To study the effect of a disordered edge, 
lattice sites from both outer edges are removed [9] [33] 
with a probability p e = 0.3 and an average over twenty 
edge disorder realizations was carried out. The average 
conductance for a nanoribbon of width W = 15.9 nm, 
d e = 8 nm and L = 50 nm is shown as the continuous 
line of Fig. 6(a). Comparing with the dashed line (same 
nanoribbon without outer edge disorder) we find that the 
anti-resonances of the first channel are robust structures, 
while the peaks higher channels are completely washed 
out, except for those three peaks in the second electron 
plateau. This robustness is associated, as already shown, 
by the fact that they have a high current density along 
the defect. In Fig. 6(b) we show, as a continuous line, the 
average conductance of the same nanoribbon, but moving 
the defect toward one edge, d e = 5.8 nm, comparing with 
the dashed line (same nanoribbon without edge disorder) . 
The main features are completely washed out except for 
the peaks which have a high current density along the 
defect. These results suggest that several conductance 
fingerprints, introduced by the linear defect, are robust 
against outer edge disorder showing that these extended 
defect structures could be used in real mesoscopic devices 
as a quantum wire. Such devices could be an interest- 
ing alternative to other potential modulation proposals 
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V. FINAL COMMENTS 

In summary, we calculate the transport properties of 
a linear defect, recently observed in bulk graphene |16j . 
embedded in a zigzag nanoribbon. Our model calcula- 
tion, when compared to first principle results [16] , reveal 
that the main qualitative features have been captured. 
Since the bond lengths around the linear defect are mod- 
ified, a further consistency check has been undertaken, 
namely modifying the hopping parameter t according to 
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the change in the bond lengths [35]. Such procedure re- 
vealed no significant changes in our results [36]. Look- 
ing therefore to the conductance, a variety of interest- 
ing fingerprints emerge. One of them is the presence of 
anti-resonances for electrons and bona fide Fabry-Perot 
oscillations for holes. Such features should be still diffi- 
cult to observe considering present device quality. More 
robust fingerprints, having in mind possible experimen- 
tal observation should be three fold: (i) presence of an 
even conductance plateau, (ii) a wide anti resonance in 
the first hole like plateau and, (iii) new positions of the 
conductance plateaus onsets, for the redefinition of the 
effective nanoribbon width by the linear defect, justify- 
ing the title of the present paper, that defines the lin- 
ear defect as third edge for the nanoribbon. Indeed, the 
analysis of current density distribution and conductance 
calculations in presence of outer edge disorder suggest 
that many conductance fingerprints, related to appre- 
ciable current density around the linear defects, are ro- 
bust against disorder which would be expected in real 
samples. Such samples could be obtained by means of 
ion beam irradiation [18] to create the extended defects, 
while the nanoribbon structure around the defect would 
be achieved by atomic force microscopy nanolitography 
0. 
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